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Abstract 

The protonic conductivity in the hydrogen bonded chains is investigated 

theoretically in the framework of the two-stage transport model. The strong 

interactions with optical phonon stretching mode are considered. We obtain 

the transition from the insulator to the metal-type state from the temperature 

dependencies of the hopping conductivity and analyze the influence of the 

specific Grotthuss mechanism on the transition character. We investigate 

also the main peculiarities in the frequency dependencies of the band and 

hopping conductivity parts which appear due to above-mentioned two-stage 

process of the proton migration along the chain. 
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I. INTRODUCTION 



The investigations of transport phenomena in hydrogen-bonded materials attract much 
attention due to their fundamental role in a variety of systems which are important for 
industrial, biological and other applications. It is well known that electric currents in these 
compounds are carried out almost entirely by mobile protons. On the one hand, this refers 
to proteins and biomembranes with one-dimensional channels for proton transport (for in- 
stance, proton pathways in bacteriorhodopsin) in which the proton transfer is directly in- 
volved in many biochemical processes^. Such a materials are well known as sufficiently 
well protonic conductors with the conductivity in chain direction which exceeds 10 3 times 
the conductivity in perpendicular direction!. On the other hand, superprotonic phases dis- 
covered in some H-bonded systems are characterized by drastic increase of the protonic 
conductivity up to the values about 0.1 Q ■ cm . It is generally accepte that the 
two-stage conduction mechanism is required to sustain the proton transport. The intrabond 
proton tunnelling along the hydrogen bridge is connected with the transfer of ionic positive 
and negative charged defects, whereas the intermolecular proton transfer due to the reori- 
entations of molecular group with proton leads to the breaking of the hydrogen bond and 
creation of a new one between another pair of molecular complexes (so-called orientational 
Bjerrum defects). It should be noted that the formation of the hydrogen bridge induces the 
distortion of groups involved in hydrogen bonding towards the proton that results in the 
shortening of this bondi. By this means the protonic polaron is localized between distorted 
ionic groups at sufficiently low temperatures, giving rise in particular to the dimerized low- 
temperature phases of some superprotonic conductors^. As has been shown in Ref. || the 
small polaron is formed due to the strong coupling of proton with the optical stretching 
vibration modes of the molecular groups. 

Since the transport process has a two-stage character, the problem of the analysis of each 
stage contribution to the total proton mobility is of particular interest as far as there exist a 
systems with different relative concentration of ionic and orientational defects. It should be 
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pointed out that whereas in systems like KH2PO4 a rapid dielectric relaxation is attributed 
to an abundance of mobile ionic defects and the rotation is sluggish, in ice the orientational 
D and L faults dominate the dielectric relaxation due to their greater numbers that limits 
the conductivity in ice. Indeed, the latter fact should manifest itself in the peculiarities of 
measured conductivity, and careful theoretical investigation could provide a possibility to 
support existing conjectures concerning the transport mechanism. 

We start in this work with a double-minimum model for the protonic potential, which is 
modified by strong interaction with lattice vibrations. Aside from the proton motion in a 
double-well potential within the bond, we assume the interbond proton hopping which makes 
possible migration of protonic polaron along the chain. We neglect here the interproton 
correlation effects that can be valid in the case of low proton concentration in chain. However, 
we take into account the possibility of proton exchange between our selected chain and 
surrounding. 

We study in this work the influence of ionic group displacements on the proton subsystem 
behavior. It is worthy to note that since the proton contribution into the electric current 
dominates in considered here H-bonded materials, the interaction of proton subsystem with 
the lattice distortions may lead to localization of protons and the system would therefore 
be an insulator. However, as we show in this work, the temperature dependence of the 
conductivity a changes from that typical of an insulator (decrease of o upon lowering T) to 
that typical of a conductor (increase of a upon lowering T). This effect points out to the 
existence of the so-called "metal" -insulator transition which was observed in the electronic 
semiconductors!. Moreover, the proposed here theoretical approach enables us to analyze 
the influence of the specific transport mechanism on the transition character. We show that 
in our case the transition to the " metallic" state with temperature increase has the peculiar 
two-step character which attributed to the two-stage nature of the proton migration pro- 
cess. Besides the specific temperature dependence, the above-mentioned complex transport 
mechanism manifests itself in the frequency dependencies of a. Although we neglect initially 
the correlations between protons, it appears from the dynamical conductivity behavior, that 
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the proton motion in the chain is a strongly correlated process (the neighbouring protons 
tend to move in pairs along the chain). While the first step of our analysis consists in the 
quasi-one-dimensional chains study, we believe our results can also be relevant for other 
hydrogen-bonded materials. 



II. DESCRIPTION OF THE MODEL 

The object of our consideration is the infinite chain shown in Fig. [l|(a). However, to avoid 
the geometric complexities introduced by the kinks in such a zig-zag chain, we consider in 
our model linear chain (see Fig. |l](b) where two neighboring chains are shown). The process 
of the proton transfer in the double-well H-bond potential is represented as the quantum 
tunnelling between two proton states with the intrabond transfer integral Qt 

tt T J2( C taClb + ct b Cla), (1) 
I 

where c^, c\ v denote proton creation and annihilation operators in the position (I, v = a, b) 
of the chain. Besides that, we describe the interbond reorient ational proton hopping in 
two-level approximation as the quantum tunnelling effect with the hopping amplitude Qr 

flfl£(ci!; lia c { 6 + 4q J+ i >a ). (2) 
l 

In this way within the framework of orientational-tunnelling model proposed in Ref. [K] the 
two-stage proton migration mechanism can be considered as the sequential migration of the 
ionic and orientational defects. The proton energy spectrum and dielectric properties of the 
finite H-bonded complexes were analyzed on the basis of exact diagonalization procedure^. 

As far as such a double-well chain is just a structural component of the system we also 
admit a possibility of proton exchange between the chain and surroundings by considering 
the system thermodynamics in the framework of the grand canonical ensemble with inclusion 
of the proton chemical potential 

-H^niv (3) 

l,v 



which is to be determined at the given proton concentration in the chain from corresponding 
equation for the chemical potential. 

Our main interest here is to analyze the influence of the longitudinal optical ionic group 
vibration modes on the proton subsystem transport properties. We consider the optical 
in-phase stretching vibrations of ionic groups in the chain which induce their displacements 
with respect to the surrounding chains as identified in Fig. |l|(b) by dashed arrows. The 
interactions of protons with this optical mode cause the difference of the potential minima 
(I, a) and (I, b) depth within the bond 

Y. T M)i n ia-ni^ b ){b q + bt q ). (4) 

l,q 

Here 77(g) = gJh/2MNuj(q) exp[iqld] where g is the corresponding coupling constant, M is 
the effective ionic group mass, N denotes the number of hydrogen bonds in chain and d is 
the lattice parameter. The optical phonon branch creation and annihilation operators are 
denoted by b't and b q respectively. Furthermore we assume the harmonic approximation for 
the lattice vibration energy 

X>fa)tfV (5) 

<i 

First of all let us consider the case of the isolated chain without a coupling to the phonon 
bath. Since the Hamiltonian ([]])- © can be exactly diagonalized, the proton energy spectrum 

£„(*) = ±|t fc |, t k = n T + n R e-* kd (6) 

forms two energy bands with the bandwidth Ae = Qt + — \ &>t — Qr \ ■ The energy gap 
in this case is = 2\Q<r — Eliminating one of the elementary transport process 

by setting the hopping amplitude Qt = or Qr = 0, we can see that both of the energy 
bands degenerate into the two energy levels and the quantum fluctuations between these 
two system states could be derived. It is clear that in the case when n = 4p J2( n iv) = 1 (one 

lv 

proton is averaged within the bond) the lower band is filled and the chemical potential /x is 
centered between bands - thus the material is in insulator state. However, for n — |, as an 
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example, only half of the lower band is filled and this corresponds to the case of protonic 
conductor that occurs for example in superionic phases of superprotonic crystals. 

We will discuss afterwards the consequences of the proton-phonon coupling effect focusing 
on the analysis of the strong coupling regime. Since the low-order perturbation theory is 
insufficient in this case, we apply a canonical Lang-Firsov transformation 

U = e* s , S = "£n lv v lv , Ujo = - Uj _ li6 = i^^-(6,-6+). (7) 

w V M?) 

We present the transformed Hamiltonian 

H = e lS He~ iS = H + H t , (8) 

where 



Hq = -ft n i" + E o H ni a ni-i,b + hv(Q)bqb q , (9) 
lv l q 

I I 

Here ft — /i + Eq is the proton chemical potential which is renormalized due to the lattice 
polarization and the formation of the protonic polaron as the result. The notation E = 

J2 \ri(q)\ 2 /hu(q) = {Tig) 2 /2M(hu ) 2 denotes the binding energy of the small polaron which 

q 

forms in the strong proton-phonon coupling regime and the phonon operators 

_ ynM+r l+1 ( q ) + _ 2 y^(6„-fe+) 
& T = e " , $ l R = e " , (10) 

give rise to the protonic polaron band narrowing factors with their thermal averages given 
by 

-£ll(9)l 2cth 3/JM9)/(M9)) a -2^|r ; (g)| 2 cthl/3^(g)/(^(g)) 2 

(<J>T> = e « ■ , (& R )=e « ■ =(^) 2 . (11) 

Furthermore we consider the case of n — 1 with the filled lower proton band which 
corresponds to the insulator state of the system. 
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III. CURRENT OPERATOR AND OPTICAL CONDUCTIVITY 



To analyze the protonic conductivity along the chain direction, we should evaluate the 

following expression in the framework of the linear response Kubo theoryEl 

1 roc r/3 

a(u,T) = — dtexp\i(u + ie)t] dX(J(t-ihX)J(0)), (12) 
Nd Jo Jo 

where the proton current operator J = [H, x] = Jt + Jr [x = n w r iu is the proton Polar- 
is 

ization operator and r\ v is the proton coordinate in the z/th position of the /th bond). After 
the transformation (|7|) the parts Jt and Jr which represent the intrabond and interbond 
proton transfer respectively, have the following form 

Jt = ^R ab J2(& T cic lb - ^c+c la ), (13) 

JR = ^RrJ2(&Rct hb C la - ^C+Q_ 1)6 ), 

where R a b is the separation between the potential wells within the H-bond and R r is the 
distance between the nearest proton positions on the neighbouring bonds adjacent to the 
ionic group. In this case the polaronic bandwidths As = Qt+^r— \&t~ Rr| (^t — ^t(^t)> 

= Qr(& r )) narrow exponentially with temperature increase due to the proton-phonon 
scattering processes. Because of this, the polaronic bands degenerate practically to localized 
states and the proton migration proceeds as the hopping of the protonic polaron between 
localized positions in the lattice at high temperatures. 

It is well known that the total polaronic conductivity can be represented as a sum of 
two parts0 : a = a co h + o hov . The band conductivity a coh dominates at low temperatures, 
whereas the hopping term ah op makes the main contribution at high temperatures, when the 
bandwidths are small and the particles become trapped at individual sites due to phonon 
interactions. We consider first the protonic low-temperature band conductivity. Since the 
scattering caused by phonon interactions can be neglected in this case, the proton transport 
is assumed to be the motion in the thermally averaged phonon field. Then evaluation of the 
correlation functions in (O) using the Wick's theorem yields the following expression for 



the real part of a co h 



cr C oh = cr° coh + af oh , (14) 



where the Drude-type term can be written as 

e 2 (3 ^(n k tt-n* k t k ) 2 



with Q k = VtrRab - fi R i? r e ikd , t k = Q T + Q R e ikd , and af oh is given by 

^ ab _ e 2 1 [(n 2 T - n 2 R )(R ab + R r ) + {hu/2)\R ab - R r )f 

Gcoh fid (hu) 2 ch/3/2 + ch[3huj/2 JqCItCIr) 2 - [(hcu/2) 2 - (fl T + f^) 2 )] 2 

Here the notations (n ku ) = [1 + exp(/3(±|i&| — /i))]" 1 denote thermally averaged Fermi 
distributed polaronic band occupancies. We see that the first term a® oh corresponds to the 
intraband transfer processes giving rise to the absorption peak at uj = 0, whereas a^ h is 
associated with the interband transfer processes and leads to appearance of the van-Hove- 
type singularities in the absorption at %uj = ±(f2 T + Q R ) and Tiuj = ±|f2y — Qr\. 



IV. THE PHONON- ACTIVATED PROTON TRANSPORT 

Furthermore, we will focus on the analysis of the polaronic transport which corresponds 
to phonon-activated hopping, because of the dominant contribution of this part to the total 
protonic conductivity at high temperatures. 

We consider the following correlator in the conductivity relation (O) which characterizes 
the hopping transport 

K hop = (J(z)J(0)) hop = K T hop + Kl p + K™ (17) 

where z = t — ihX and 

n IV 

(l>^(z) C +(z)q_ 1)6 (z)$^+_ li6 c, a )) , (18) 
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n IV 

(¥*(z)cl(z)c la (z)¥ R 4 a c v ^ b ) + 1;b (z)c la {z)^c+ a C Vb ) + 

Here the operators = — (& a ) (en = T, R) correspond to the phonon scattering processes 
which accompany the proton transfer. We perform evaluation of using the Hamiltonian 
with averaged values of the transfer integrals Q a . This allows us to decouple the proton- 
phonon correlation functions 

(KizKizMz^Z 1 ^^) - (^{z^H^zMz)^^) (19) 

into proton and phonon parts which can be evaluated separately. The processes that we 
are interested in are going on at considerably high temperatures and in the strong proton- 
phonon coupling regime (we will see below that the strong polaron effect actually occurs in 
these systems). This allows us to assume 

V »i (20) 

Y {huj{q)) 2 shf3hu{q)/2 1 ' 

We take a dispersionless approximation for the phonon frequency: uj(q) = uq. In this case 
E = (kg) 2 /2M(huJo) 2 and evaluation of the phonon correlation functions yields 

^ RR = (& R (z)& R *) = (^(z)& R ) = e-^e-^ h ^ 2 / Af n5w, 

e -^cth^o/2 e -(, + ^/2) 2 /4f| (fe + j ) 



The parameters introduced in (|2T| ) 



tI = fi = — — , r' 2 = — — (22) 
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characterize the average hopping time lengths between two localized positions in the chain 
(the hopping time within H-bond, the interbond reorientational hopping time and the time 
for pair polaron transfer within the neighbouring bonds respectively). 

Next we consider the proton part of the correlation functions. We introduce the k- 
representation for the Fermi operators 

and further perform the following unitary transformation 

Cka = J- {[tk ~ i\t k \]Cka ~ + itkjCkb) , 

®>k 

Ckb = -j- ([\tk\ - it* k ]c ka + [t* k + i\ik\]c k b 
dk v 



where dk = 2y \tk\y \tk\ — ^st k . Then the polaron Hamiltonian can be reduced to the diagonal 
form Y.ku^kunkv {e kv = -p, ± \t k \) and the correlation functions {cX v {z)c k 'v\z)cX xVx c k ^ can 
be evaluated using the Wick's theorem (in this case we use the Schrodinger representation 
for the polaron operators c kv [z) = e~ %£kvZ / h c ku ). 

We use procedure proposed ini with the deformation of the integration contour in the 
complex plane for the integration of expressions flTSP over t and A. Let us present the 
resulting expression for the real part of the conductivity that will be analyzed afterwards 

a hop (uj,T) = cr hop (u,T) + Y,c r hop{u,T), (23) 

vv' 

n° d ) t\ — 2e2 V^ sh ^/ 2 v Ur P - f > 2 -i- r> p--T- 2 wrf*f ^ 

VhopK^,! ) - , N2h 2 n /9 Z^U G ^ e +C T e )M(t k t k ') + 



d.XVi 1 huj/2 

\tk\\tk' 

t) = g {^"^ + ^ e ~^ 2 cos ^ - k '^ d+ 

C R e~^ 2 + C TR e- f T" 2 (cos kd + cos k'd)} e /3(£fc ^ £ ^' )/2 (n^)(l - (n feV )). 
Here cD = w + £k »~£ k '»' and 



C T = ^^e-^/ 2 , C' T = 2n 2 T R 2 ab f' T e-^e 
C R = n R R 2 f R e~^\ C TR = -2n T n R R ab R r f T e-^/ 2 e~^ cth ^ 0/2 
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The first term in (|23|) appears in the result of pairing of the operators c£ v with Ck' V ' in the 
proton correlation functions with equal values of time argument, whereas another terms a^ p 
arise due to the pairing of such operators with different time arguments. 

Let us discuss first the frequency dependence of the hopping conductivity. From expres- 
sion fl23] ) it follows that a^op has several peaks which point to the important contribution 
of various elementary dynamical processes existing in the system. It is possible to estimate 
analytically in the limit (3Tiu » 1 (in the high frequency range or at lower temperatures) 
the characteristic frequency values at which the conductivity maxima exist. The first maxi- 
mum at u>i ~ ^pj corresponds to the pair phonon-activated hopping within the neighbouring 
H-bonds. Another maximum at co>2 ~ Jpr arises due to the intrabond polaron dynamics, 
whereas the peak at u 3 ~ is associated with the reorient at ional hopping of the protonic 
polaron between bonds. The frequency dependencies of the real and imagine parts of the 
hopping conductivity are shown in Figs. 0, 3 for different temperatures and values of polaron 
binding energy E . We observe three distinct peaks at sufficiently low temperatures, which 
broaden significantly when temperature increase. Finally, at a considerably high temper- 
atures the behavior <Jho P {u) has the Drude-type character with maximum at u = and 
smoothing decrease with increasing frequency. 

The two-stage transport mechanism manifests itself in the two peaks at u = U2 and 
uj = 003 observed in the frequency dependence of the optical conductivity which shift to 
the high-frequency region when the polaron binding energy E increases (see Fig. 3). With 
the increase of E the conductivity lowers and the maxima broaden due to the stronger 
localization of protonic polaron in one of the potential well of the hydrogen bond. We note 
that the values of the hopping integrals Qt and CIr practically do not influence on the 
positioning of these maxima but affect strongly the redistribution of the maxima heights. 

Figure 4 shows the comparison of theoretical conductivity obtained for lithium hy- 
drazinium sulphate Li(N 2 H 5 )S04 (LHS) using (|23| ) with experimentally measured values@ 
at the frequency 500 KHz. This compound consists of the chains of hydrogen bonded hy- 
drazinium ions running along the [001] direction and can be considered as the example of 
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quasi-one-dimensional protonic conductor^. In this case we obtain hooo = 800 cm -1 and 
E Q = 0.98 eV. According to our results, the dominant contribution to the conductivity 
in this frequency range arises from the intrabond hopping with Qt = 0.04 eV. The rela- 
tion Eq > 2Qt follows immediately from the values of our parameters and is valid for the 
small-polaron regimeS. 

Let us analyze the theoretical hopping conductivity in a wide temperature interval for 
the various system parameters. Since the proton is localized in the chain due to the strong 
polaronic effect, we expect intuitively that change in an external control parameters may 
have a determining influence on the proton state and lead to its derealization. In particular, 
it is of prime importance to analyze in more details the influence of temperature. On the one 
hand, it is well known that the effect of temperature is dramatic in the systems which undergo 
the transitions to superprotonic phases, which accompanied by significantly increase of the 
protonic conductivity. In this case the proton state changes especially drastically: in the low- 
temperature ordered phases the protons are localized in well-defined positions of the lattice 
and can migrate with sufficiently high activation energy (~ 1 eV) due to the strong coupling 
with displacements of ionic groups participating in hydrogen bond formation, whereas in 
the superionic phases the transport phenomenon is connected with dynamically disorder of 
the hydrogen bond network, resulting in an increase of the number of possible positions 
for protons which can hop between them with much lower activation energy (~ 0.1 eV). 
On the other hand, the influence of temperature on the localized state of protonic polaron 
can be especially important in a variety of protonic conductors, both in inorganic systems 
and biological membranes. Because of this, we consider the temperature dependence of the 
hopping conductivity at the given frequency of external field. 

The real part of the Oh ov as a function of temperature has peculiarities which arise due 
to derealization effects. We consider several characteristic frequency regions: u ~ uj 2 (here 
the intrabond hopping contributes mainly to the total conductivity), u ~ U3 (when the re- 
orientational interbond hopping dominates) and uo = (static dc conductivity). The effect 
of two-stage transport process shows itself most clearly in the high-frequency region. The 
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first case is shown in Fig. 5. The two maxima are observed which shift slightly towards the 
higher temperatures when E increases. It is easy to understand the origin of this nonmono- 
tonic behavior of Ohop from Fig. 6, where the temperature dependence of a hap is compared 
with two extreme cases that one of the transport component (interbond or intrabond) is 
eliminated. Now we can conclude immediately that the first sharp increase of a^op with 
further rapid lowering originates from the intrabond hopping, whereas the second smoother 
maximum on higher temperature interval corresponds to the interbond contribution. The 
raising of temperature suppresses the polaron localization effect: the thermal proton dynam- 
ics leads to the detrapping of protonic polaron which accompanied by drastically increase 
of crhop- However, with the further temperature increasing the conductivity lowers again 
due to the multiphonon scattering processes suggesting that the system undergoes a transi- 
tion to the metallic-type state. Therefore, we can conclude that in our system the possible 
transition from the insulator state has a two-step character: first the protonic polaron is 
delocalized within the H-bond, and only with further temperature increase it becomes de- 
localized between bonds. The situation changes alternatingly in the vicinity of lo ~ (see 
Figs. 7 and 8). We see that in this case first the tendency for interbond derealization oc- 
curs, whereas the intrabond contribution dominates at higher temperatures. The last case 
&ho P (u = 0) is shown in Fig. 9. for different relations between Qt and Qr. We see the 
single maximum which corresponded to the transition to metallic-type state appeared due 
to the dominant contribution of reorientations (when Qr > Qt, Fig. 9(a)) or intrabond 
hopping (when Q T > Q R , Fig. 9(b)). In the case Qt ~ (Fig. 9(c)) the both parts con- 
tribute practically equally to the transition occurrence. It is worth noting that despite of 
the single smooth peak, we can also distinguish two steps in the dc conductivity behavior 
which becomes clear after comparison of our curves with the two extreme cases (Qr = and 
Qt — 0). We can conclude that the protonic polarons become delocalized between bonds 
at lower temperatures, and the intrabond derealization proceeds on further heating. This 
effect can be explained by the fact that the reorient at ional hopping has higher activation 
energy than the interbond motion (see expressions (|23|) for the conductivity) and thus the 
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first type of hopping has the tendency to damp at lower temperatures due to the multi- 
phonon scattering than the second type of hopping. Thus our model allows us a possibility 
to analyze thoroughly the effect of either of two transport stages on the total conductivity 
picture that can be verified by experimental measurements. 

V. CONCLUSIONS 

In the present work the protonic conductivity of the quantum quasi-one-dimensional 
hydrogen-bonded chain has been studied in the framework of the two-stage orient at ional- 
tunneling model. The interaction of protons with the optical stretching mode is considered. 
The expression for the conductivity coefficient as a function of temperature and frequency 
is obtained using the Kubo formula for various values of hopping parameters and polaron 
binding energy. The evaluated hopping conductivity is in good agreement with the ex- 
perimentally measured values for LHS system. This allows us to conclude that the strong 
proton-phonon coupling regime exists in this type of systems. The analysis of theoretically 
obtained temperature dependencies of Oh op in wide temperature interval revealed that the 
transition from the insulator state (increase of a with increase of T) to the metal-type state 
(decrease of a with increase of T) occurs in system on heating. The specific two-stage Grot- 
thuss transport mechanism manifests itself in the two-step character of the transition as 
well as in the two distinct peaks of dynamical conductivity. The obtained theoretical results 
allow us to believe that the proposed theory should also apply to other hydrogen bonded 
systems and to investigate the consequences of assumptions of the Grotthuss mechanism for 
proton migration in hydrogen bonded compounds. 
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FIGURE CAPTIONS 



Figure 1. (a) Zig-zag hydrogen-bonded chain, arrows indicate the possible path of proton 
migration along the chain, (b) Simplified model chains, the in-phase displacements of ionic 
groups identified by the dashed arrows. 

Figure 2. Frequency dependencies of the hopping conductivity c%, p (u;) = o hov jc^ 
(co = u = Lo/uo) for E = E /huj = 2 and different temperatures r = kT/hu> ; 

R a b/ r o — 1-3, R r / r o = 1-6. ojt = (It — 0.8 and uj r = Q R /hoj = 0.5 are reduced 
hopping amplitudes. The real and imagine parts are indicated by bold and thin curves 
respectively. 

Figure 3. Frequency dependencies of 'Sldhop at r = 0.25 for different values of E Q and 
hopping amplitudes. 

Figure 4. Comparison of the temperature dependencies of the protonic conductivity, mea- 
sured for the crystal LHS and evaluated using (|23|) . 

Figure 5. Temperature dependencies of tft&hop for different values of E at Co = 5.5, ujt = 0.2 
and ujr = 0.6. 

Figure 6. Comparison of 3t,ah op (T) for E = 2.5, uj = 5.5, u>t = 0.2 and ujr = 0.6 with two 
extreme cases when the intrabond or interbond contribution is eliminated. 
Figure 7. Temperature dependencies of $t.ahop for different values of E at uj = 8.7, ujt = 0.6 
and ujr = 0.3. 

Figure 8. Comparison of $lahop(T) for E = 2.3, uj = 8.7, ujt = 0.6 and uj r = 0.3 with two 
extreme cases when the intrabond or interbond contribution is eliminated. 
Figure 9. Temperature dependencies of dc hopping conductivity for different values of ujt 
and ujr at Eq = 2.5. 
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